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The dynamics of magnetosphere-ionosphere coupling has been investigated by means of a 
two-dimensional two fluid MHD model including anomalous resistivity. When field-aligned 
current is generated on auroral field lines, the disturbance propagates towards the ionosphere in 
the form of a kinetic Alfven wave. When the current exceeds a critical value, microscopic tur- 
bulence is produced, which modifies the propagation of the Alfven wave. This process is 
modelled by a non-linear collision frequency, which increases with the excess of the drift velo- 
city over the critical value. Turbulence leads to absorption and reflection of the Alfven wave, 
partially decoupling the generator from the ionospho.,:. The approach to a steady-state is 
strongly dependent on the presence or absence of the turbulence. The current is self-limiting, 
since a current in excess of critical causes a diffusion of the magnetic field perturbatioki and a 
reduction of current. The transverse scale size is determined by a balance between the non- 
linear steepening of the wave and the diffusion caused by the turbulence. 

The coupling of magnetosphere and ionosphere by field-aligned currents is one of the 
most important problems of magnetospheric and auroral physics. Previous models of auroral 
electrodynamics generally either have ignored microscopic processes such as plasma turbulence, 
or have considered only the microscopic processes, without reference to the large-scale 
processes which feed energy to the turbulence. We have attempted to bridge this gap by means 
of a two-dimensional MUD model, which includes a term to represent a resistivity due to 
plasma turbulence. 
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Thu computer model allows variations in the xt plane, where x represents the latitudinal 
direction and z, the direction of the background magnetic field, which is a^umed to be con- 
stant in space and time. The Held alignwl current is, then, represented by a magnetic perturba- 
tion By. The electric field then has components in the xz plane. The equations describing 
induction, ion polarization drift, continuity, and electron motion along the field, in conjunction 
with Ampere's Law, then determine the evolution of the system. In the linear limit, these 
equations yield the dispersion relation for kinetic Alfv^n waves in a cold plasma: 

In the model, terms which are less than the electron-ion mass ratio are neglected, which leaves 
the convective nonlinearity as the main nonlinear contribution. In this limit, electron motion is 
along the background Reid, whereas the ions move only across field lines. 

Turbulence is modelled by means of an effective drag term introduced into the electron 
equation of motion. When the electron drift velocity m is greater than the critical value for ins- 
tability, this drag term h " he form: 
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The drag term is of the opposite sign for u < -Uc, and is zero when |u| < For clectrof- 
talic ion cyclotron turbulence saturated by resonance broadening, one finds « 2, vi =« 0.2n,, 

vj =« 0.04n,, where 0, is the ion gyrofrcqucncy. In the model, vj and vj .tre taken to be 
independent of space and time, while is allowed to be a function of ;. 


On one end of the system, representing the ionosphere, a height integrated Pedersen con- 
ductivity is specified which causes Alfvdn waves to be reflected. For realistic ionospheric values 
(I^ — 10 mhos) the Alfven wave is almost perfectly reflected. On the other end, a generating 
field is imposed. Both current and voltage generators can be considered by specifying Dy and 
respectively, as a function of .x and of time at this end of the system. 


In the absence of turbulent drag, the system evolves much as in the work of Cmv/z and 
Boswell (1979). Since the ionospheric height integrated conductivity is generally much greater 



3 


than the Alfvin wave impedence - c*/ (4ir + klU an Alfvdn wave pul» 
incident on the iont^phere will be reflected with Eg reversed and By in the same direction. If 
Eg is fixed in the gencictor, the current will increase until a steady state is achieved in which Eg 
is indei^ndent of z and the ionospheric Pedersen current is given by the imposed Eg and the 
assumed Peder^n conductivity. Similarly, if By is imposed. Eg will decrease until Ohm’s Law 
is satisfied in the ionosphere. 

In the presence of turbulence, the situation is changed dramatically. When an Alfvdn 
wave pulse encounters a region of turbulence, it is partially reflected, transmitted, and ateorbed 
[Lysak and Carlson, 1981]. The ;vave pulse is reflected in the opposite sense a.s reflections from 
the ionosphere; that is, the By comi>onent is reversed. Thus regions above a turbulent region 
will have strong Eg and weak By components, while on quiescent field lines, Eg will be weak 
and By strong. In this case, if Eg is imposed, magnetic field lines will diffuse on field lines 
where turbulence is present, leading to a broader distribution of current than that originally 
imposed. On the other hand, when By (and thus, the current profile) is imposed. Eg is 
enhanced in regions where the current, and thus the turbulence, is strongest. This leads to the 
formation of narrow V-shaped potential structures on field lines where turbulence is present. 

The effect of ionospheric and turbulent regions on the generator may be discussed by 
recalling that Eg may be interpreted as a plasma flow in the y direction. When a pulse is 
reflected by the ionosphere, the Acid-aligned current is increased, causing an increase in the 
j X B force which resists the plasma flow. When the pulse reflects off the turbulent region, 
however, the generator is cut off from tne ionosphere and the flow speed may increase due to a 
reduction in the j .< B force. This may provide a mechanism for plasma injection on field lines 
which contain turbulence. 

The approach to a quasi-steady-statc can occur much more rapidly in the presence of tur- 
bulence. V-shaped potential structures are formed within one or two bounce periods. Integral 
properties of the system, such as potential drop and total energy dissipation, may reach constant 
values while local properties are still oscillating. Thus, the presence of a steady arc structure 


does not necessarily imply that the flux tube is locally in a steady state. 
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